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ABSTRACT 

Simulations in general relativity show that the outcome of collapse of a marginally unstable, 
uniformly rotating star spinning at the mass-shedding limit depends critically on the equation of 
state. For a very stiff equation of state, which is likely to characterize a neutron star, essentially 
all of the mass and angular momentum of the progenitor are swallowed by the Kerr black hole 
formed during the collapse, leaving nearly no residual gas to form a disk. For a soft equation of 
state with an adiabatic index F — 4/3 <C 1, which characterizes a very massive or supermassive 
star supported predominantly by thermal radiation pressure, as much as 10% of the mass of the 
progenitor avoids capture and goes into a disk about the central hole. We present a semi-analytic 
calculation that corroborates these numerical findings and shows how the final outcome of such 
a collapse may be determined from simple physical considerations. In particular, we employ a 
simple energy variational principle with an approximate, post-Newtonian energy functional to 
determine the structure of a uniformly rotating, polytropic star at the onset of collapse as a 
function of polytropic index n, where F = 1 -I- 1/n. We then use this data to calculate the mass 
and spin of the final black hole and ambient disk. We show that the fraction of the total mass 
that remains in the disk falls off sharply as 3 — n (equivalently, F — 4/3) increases. 

Subject headings: black hole physics - relativity - hydrodynamics ~ stars: rotation 



1. Introduction 

Determining the final state of a rotating star 
undergoing gravitational collapse is an important 
issue in relativistic astrophysics. Such a collapse is 
the principle route by which a rotating black hole 
forms in nature. Several different scenarios involv- 
ing the coUape of rotating stars to black holes have 
been the focus of considerable attention recently. 
The collapse of massive stars in hypernovae ex- 
plosions ( "coUapsars" ) may be the origin of long- 
period gamma-ray bursts (MacFadyen & Woosley 
1999; MacFadyen, Woosley & Heger 2001). The 
remnant of a binary neutron star, following in- 
spiral and coalescence, is likely to be rapidly ro- 
tating and undergo collapse, either prompt or de- 
layed, depending on its mass (Shibata & Uryu 



2000, 2002; Baumgarte, Shapiro & Shibata 2000). 
Binary neutron star merger and collapse provides 
a plausible scenario for generating short-period 
gamma-ray bursts (Narayan, Paczynski & Piran 
1992; Ruffert & Janka 1999) and a good candi- 
date for the detection of gravitational waves by 
Advanced LIGO and other high-frequency, gravi- 
tational wave laser interferometers. The collapse 
of a rotating supermassive star may be a promis- 
ing route for forming the seeds of supermassive 
black holes (see Shapiro 2003 for a review and ref- 
erences), which are likely to reside at the centers 
of many, and perhaps most, bulge galaxies (Rich- 
stone et al. 1998; Ho 1999), including the Milky 
Way (Genzel et al. 1997; Ghez et al. 2000; Schodel 
et al. 2002) and are believed to be the engines that 
power active galalactic nuclei (AGNs) and quasars 
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(Rees 1998, 2001). 

General relativity induces a radial instability 
in compact stars to catastrophic collapse. Rota- 
tion can support larger masses in stable equilib- 
rium, but all stars become unstable to collapse 
when they are sufficiently compact. Uniformly ro- 
tating stars can have masses exceeding the maxi- 
mum allowed mass of a nonrotating spherical star 
by about < 20%. Cook, Shapiro & Teukolsky 
(1994a,b) have constructed numerical solutions of 
such "supramassive" stars, both for polytropes 
and for 14 realistic candidate nuclear equations 
of state. Differentially rotating stars can be "hy- 
permassive" , with masses exceeding the maximimi 
value for supramassive stars by factors or two 
or more (Baumgarte, Shapiro & Shibata 2000; 
Lyford, Baumgarte & Shapiro 2003). Viscosity 
(molecular or turbulent) and magnetic fields drive 
differentially rotating configurations to uniform 
rotation on secular timescales. [See, e.g., Shapiro 
2000 and Cook, Shapiro & Stephens 2003 for re- 
cent Newtonian calculations of magnetic braking 
and viscous damping of differential rotation, and 
Liu & Shapiro 2004 for approximate, general rela- 
tivistic calculations of these phenomena; see Duez 
et al. 2004 for detailed, general relativistic, numer- 
ical simulations of viscous damping of differential 
rotation.] 

Here we focus on uniformly rotating stars and 
consider the collapse of supramassive stars. The 
first investigation of the gravitational collapse of 
supramassive stars was performed by Shibata, 
Baumgarte and Shapiro (2000), who employed a 
fully relativistic hydrodynamics code in three spa- 
tial dimensions plus time. They studied the col- 
lapse of a marginally unstable, relativistic poly- 
trope with polytropic index n = I spinning at 
the mass-shedding limit. The mass shedding limit 
is the maximal spin rate for a uniformly rotating 
star; here matter on the equator moves in a Ke- 
plerian (geodesic) orbit about the star, supported 
against gravity entirely by centrifugal forces. Stars 
with n = 1 have stiff equations of state, are not 
very centrally condensed and are fairly relativistic 
at the onset of collapse, with Re/M « 5.6, where 
Re is the equatorial (circumferential) radius and 
M is the total mass-energy. Such stars provide 
crude models of rapidly rotating, massive neutron 
stars. The outcome of collapse is a Kerr black 
hole containing all the rest-mass Mq of the ini- 



tial configuration, and essentially all of the total 
mass-energy M and angular momentum J of the 
initial configuration (apart from the small amount 
radiated away by gravitational waves). The rea- 
son that no matter remains outside the hole is that 
at the equator the specific angular momentum of 
the gas, which is strictly conserved during axisym- 
metric collapse, is smaller than jisco , the specific 
angular momentum of a particle at the innermost 
stable circular orbit (ISCO) about the final hole. 
Hence all the interior matter in the star is captured 
by the hole. This situation is similar to the col- 
lapse of a nonrotating spherical star, which forms 
a Schwarzschild black hole without any ambient 
disk. By contrast, Shibata & Shapiro (2002) per- 
formed a fully relativistic hydrodynamics simula- 
tion in axisymmetry of the collapse of a marginally 
unstable n = 3 polytrope at the mass-shedding 
limit. Such a configuration models the likely end- 
point state of an evolved, radiation-dominated, 
very massive or supermassive star at the onset of 
collapse (Baumgarte & Shapiro 1999). By contrast 
with an n = 1 polytrope, an n = 3 configuration 
has a soft equation of state, is very centrally con- 
densed and is nearly Newtonian {R^/M 620) at 
the onset of collapse. A rotating black hole again 
forms during the collapse, but in this case the mass 
of the hole contains only about 90% of the total 
rest-mass of the system (with a spin parameter 
Jh/Mfi 0.75). The remaining gas forms a rotat- 
ing disk about the nascent hole. Here the angular 
momentum of the infalling matter in the outer- 
most layers of the collapsing star exceeds iisco; so 
the gas in the outer regions escapes capture and 
forms a disk. This result has been corroborated by 
a semi-analytic calculation by Shapiro & Shibata 
(2002), which we generalize below for arbitrary n. 

Recently, Shibata (2003) performed systematic 
numerical simulations in axisymmetry to investi- 
gate the effects of the stiffness of the equation of 
state on the outcome of the collapse of marginally 
unstable, supramassive stars spinning at the mass- 
shedding limit. He studied collapsing relativistic 
polytropes with indicies 2/3 < n < 2, represent- 
ing stiff and moderately stiff equations of state. 
He concluded that for these configurations the fi- 
nal state is a Kerr black hole and the disk mass, if 
it exists at all, is very small (< 10~^ of the initial 
mass). 

Simulations performed to date clearly show 
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that the outcome of the collapse of a maximally 
spinning, supramassive star depends critically on 
the equation of state. In this paper we employ a 
straightforward semi-analytic analysis to corrob- 
orate this conclusion and to show how the final 
outcome of these detailed simulations may be in- 
ferred from simple physical considerations. In par- 
ticular, we employ an energy variational princi- 
ple to determine the approximate structure of the 
marginally unstable, rotating star at the onset of 
collapse as a function of polytropic index n and 
use this data to estimate the mass and spin of the 
final black hole and ambient disk. We show that 
the mass of the disk depends sensitively on the 
difference 3 — n. Our semi-analytic calculations 
are based on an approximate post-Newtonian en- 
ergy functional valid for rotational energies that 
are much smaller than gravitational potential en- 
ergies. They arc most accurate for soft equations 
of state with n close to 3, since the parameters of 
the progenitor star are nearly Newtonian in this 
case and the effect of rotation in distorting the 
shape of the interior region containing most of the 
mass is minimal. Interestingly, stars with n close 
to 3 arc the most difficult to simulate with a rel- 
ativistic hydrodynamics code, since the large size 
of the progenitor star requires that any such code 
span an enormous dynamic range and integrate for 
a huge number of (Courant) timesteps in order to 
follow the collapse to completion and resolve the 
central black hole and disk reliably. In this sense, 
the two approaches - the semi-analytic analysis 
presented below and full numerical simulations - 
are complementary. 

2. The Marginally Unstable, Supramas- 
sive Progenitor Star 

In this Section we analyze the equilibrium and 
stability of a supramassive polytrope spinning at 
mass-shedding and identify the critical configura- 
tion at which radial instability sets in. Our treat- 
ment is a straightforward generalization to arbi- 
trary polytropic index n of the semi-analytic anal- 
yses of Baumgarte & Shapiro (1999). We briefly 
review the predictions of the Roche approximation 
for the rotating envelope in Section 2.1. In Sec- 
tion 2.2 we present a post- Newtonian energy vari- 
ational calculation which allows us to determine 
the critical configuration for each n and identify 
its two key parameters, Rp/M (where Rp is the 



polar radius) and J/M"^, needed to determine the 
mass and spin of the final black hole and disk. 
We calculate the mass and spin of the black hole 
and disk in Section 3. Henceforth we adopt ge- 
ometrized units and set G = c = 1. 

2.1. Review of the Newtonian Roche 
Model 

Stars with soft equations of state are extremely 
centrally condensed: they have an extended, low 
density envelope, while the bulk of the mass is 
concentrated in the core. For an n = 3 New- 
tonian polytrope, for example, the ratio between 
central density to average density is pc/ p = 54.2. 
The gravitational force in the envelope is therefore 
dominated by the massive core, and it is thus le- 
gitimate to neglect the self-gravity of the envelope. 
In the equation of hydrostatic equilibrium, 

VP 

^ = -V($ + $c), (1) 

this neglect amounts to approximating the New- 
tonian gravitational potential $ by 



In (1) we introduce the centrifugal potential 
which, for constant angular velocity f2 about the 
2-axis, can be written 

^, = -hi'{x^+y-') = -hfr^sm^e. (3) 

Integrating eq. (1) yields the Bernoulli integral 

h + ^ + $c = H, (4) 
where if is a constant of integration and 

f dP P 
h= — = {n + l)- (5) 
J P P 

is the enthalpy per unit mass. Here we have as- 
sumed a polytropic equation of state 

P = Kp^\ K, r = 1 + 1/n = const, (6) 

Evaluating eq. (4) at the pole yields 
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since ft, = on the surface of the star and $c = 
along the axis of rotation. In the following we 
assume that the polar radius Rp of a rotating star 
is always the same as in the nonrotating case. This 
assumption has been shown numerically to be very 
accurate (e.g., Papaloizou & Whelan 1973). The 
mass of the star is hardly changed from its value in 
spherical equilibrium, which to leading Newtonian 
order is given by the polytropic relation 



M = 47ri?(,3-")/(i-") 



{n + l)K 



-1 n/(n-l) 



An 



(8) 



where the Lane-Emden functions for polytropes 
appearing above are tabulated in Chandrasckhar 
(1939). A rotating star reaches mass shedding 
when the equator orbits with the Kepler frequency. 
Using eqs. (3) and (4), it is easy to show that at 
this point the ratio between equatorial and polar 
radius is 

r) -I- 

/ shedd ^ 

The corresponding maximum angular velocity is 



shedd 



3/2 



M 



1/2 



(10) 



(Zel'dovich & Novikov 1971; Shapiro & Teukolsky 

1983). 

Since the bulk of the matter is concentrated 
in the core and hardly affected by the rotation, 
the moment of inertia of inertia of the star barely 
changes with rotation and is well approximated by 
the nonrotating value 

I=^KnMRl, (11) 

where the nondimensional coefficient k„ is tabu- 
lated for many different polytropes n in Lai, Rasio 
and Shapiro (1993) and for n = 3 and n = 2.5 in 
Table 1. The ratio between the kinetic and poten- 
tial energy at mass-shedding then becomes 



T 

\W\ 



shedd 



(l/2)jn,\,,, 
(3/(5 - n)) MyRp 



= (^^j (1 - n/5) (12) 

This result predicts that T/| W| of a maximally ro- 
tating polytrope of index n is a universal constant, 
independent of mass, radius, or angular velocity. 



Inserting eqs. (2), (3), and (5)-(10) into cq. (4) 
yields the density throughout the extended enve- 
lope. 



An 



V r 
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27 R^ 



sm 



(13) 



The stellar surface is the boundary along which 
p = 0, and is thus defined by the curve r{9) given 
by 

The solution to this cubic equation is given by 
r{0) 3sin(6'/3) 



Rr> 



sin^ 



(15) 



2.2. Critical Configuration 



To determine the equilibrium and stability of 
a rotating polytrope, we express its total energy 
as the sum of the (Newtonian) internal energy U, 
the potential energy W , the rotational energy T, a 
post-Newtonian correction i?pjv and a post-post- 
Newtonian correction Eppn- Writing the terms 
in that order yields the energy functional 



E{pc; M, J) 



ksM^Pc 



(16) 



where pc is the central density and where we have 
defined j = J/M"^ and have neglected corrections 
due to deviations from sphericity. This neglect is 
justified, since these corrections scale with T/|W^|, 
which according to (12) is always very small for 
centrally condensed configurations. Even though 
the value of the post-post-Newtonian correction 
EppN is very small, this term is crucial for de- 
termining the critical, marginally stable configu- 
ration for n = 3, as emphasized by Zel'dovich and 
Novikov (1971) and Baumgarte & Shapiro (1999). 
The values of the nondimensional coefficients fcj 
are determined by quadratures over Lane-Emden 
functions. They are listed in Table 1 for n = 3 and 
n = 2.5; in our analysis below we evaluate them 
for arbitrary 2.5 < n < 3.0 by linearly interpolat- 
ing between these two values of n. 
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Table 1 

Values of the polytropic structure coefficients 



n 






k2 






k5 


3.0 
2.5 


o.issag^' 

0.27951'* 


1.4295'* 


0.63899'* 
0.67623'* 


1.204^ 
1.4202'* 


0.91829'^ 
0.86334= 


0.3312r 
0.29752= 



Note that for any polytrope, K'^/'^ has units 
of length. We can therefore introcincc nondimen- 
sional coordinates by setting K = 1 (see Cook, 
Shapiro & Teukolsky 1992). We will denote val- 
ues of nondimcnsional variables in these coordi- 
nates with a bar (for example M). Values of these 
quantities for any other value of K can be recov- 
ered easily by rescaling with an appropriate power 
of example M = K^'/^M and p = K-"p. 

Taking the first derivative of eq. (16) with re- 
spect to the central density, holding M and J con- 
stant, yields the condition for hydrostatic equilib- 
rium: 

0=— = (3/n)A;iMa;(^/"-^) - fcsM^/^ 
ox 

+ 2k:ifAf/^x - 2kiM'^/^x 
- Bk^M^x"^, (17) 

-1/3 

where x = pc ■ For stable equilibrium, the second 
derivative of eq. (16) has to be positive. A root of 
the second derivative therefore marks the onset of 
radial instability: 

= ^ = (3/n)(3/n - l)fciMa;(3/"-2) 

+ 2k3fM'^/^ - 2fc4M^/3 

- QknM^x. (18) 

To find the critical configuration at mass-shedding, 
Eqs. (17) and (18) must be solved simultaneously 
for X and M subject to the constraint 

T _ ksfM^^x^ _ kzfM^/^x 

iW^lshcdd" k,M^"^ - k, ' ^ > 

where the left hand side is given by eq. (12). 

Key parameters characterizing the critical con- 
figurations are given in columns 2 — 4 of Table 2. 



The first row for each n gives the solution found by 
solving eqs. (17) and (18) simultaneously, substi- 
tuting eq. (19). The other rows for selected n give 
the results of a careful integration of the full Ein- 
stein equations of general relativity for the critical 
equilibrium configuration. The key global param- 
eters needed to calculate the final black hole mass 
fraction and spin are the stellar compaction Rp/M 
and spin J/M"^. The value of Rp/M determined 
by the simple variational method described above 
are seen to be most reliable for nearly Newtonian 
configurations with n close to 3. Nevertheless, as n 
decreases below 3 and the equation of state stiff- 
ens, the trend toward higher compaction at the 
critical point is clearly evident. The Newtonian 
Roche model for the outermost layers yields val- 
ues for J/M^ that are consistent with the exact 
model solutions, albeit these are somewhat low, 
even for n close to 3. The discrepancy arises from 
using eq. (11) for the moment of inertia, which 
neglects the increase in / arising from oblateness 
in the outer layers, and is strictly valid only for 
extreme central mass concentration^. 

3. Calculating the Black Hole Mass and 
Spin 

In this section we generalize to arbitrary poly- 
tropic index n the iterative method of Shapiro & 
Shibata (2002) to calculate the mass and spin of 
the black hole and any disk that might form during 
the collapse of supramassive star. (They formu- 
lated the method for n = 3). Consider the implo- 
sion of matter from the envelope onto the central 

A recalculation of the exact model by Shibata (2004) for 
n = 3 yields J/M^ = 0.91, which removes much of the 
discrepancy. 
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Table 2 

supramassive star collapse vs. polytropic index n 



Critical Supramassivc Star Final Black Hole and Disk 



n 


M 


Rp/M 


J/M2 


Mh/M 


Jh/M^ 


Md/M 




3.00 


4.56 


456 


0.876 


0.89 


0.60 


0.11 






4.57'^ 


427'^ 


0.97« 


0.87 


0.71 


0.13 






4.57" 


427« 


0.97" 


0.9^= 


0.75^= 


0.1^ 




2.95 


3.83 


196 


0.602 


0.97 


0.52 


0.029 




2.90 


3.30 


119 


0.491 


0.99 


0.45 


0.011 






3.30'' 


122'' 


0.567'' 


0.99 


0.53 


0.014 




2.80 


2.54 


60.9 


0.382 


1.00 


0.37 


2.5 X 10" 


-3 


2.70 


2.02 


37.7 


0.325 


1.00 


0.32 


5.3 X 10- 


-4 


2.60 


1.65 


25.5 


0.287 


1.00 


0.29 


9.7 X 10- 


-5 


2.50 


1.38 


18.2 


0.259 


1.00 


0.26 


1.2 X 10" 


-5 




1.29'' 


28.9** 


0.389'' 


1.00 


0.39 


2.7 X 10- 


-4 



^Baumgarte & Shapiro (1999), Table 2 

''Cook, Shapiro & Teukolsky (1994a), Table 2, converting Rg to Rp 
using Rp/Re « 2/3. 

'^Relativistic hydrodynamics simulation by Shibata & Shapiro (2002) 



black hole formed from the collapse of the centrally 

condensed interior region. If the specific angular 
momentum of the imploding envelope matter is 
below iiscO) the specific angular momentum of 
a particle at the innermost stable circular orbit 
(ISCO) about the hole, the matter will be cap- 
tured. If the angular momentum of the infalling 
matter exceeds .jisco j the matter will escape cap- 
ture and continue to orbit outside the hole, form- 
ing a disk. This capture criterion is well-supported 
by the numerical simulations of rotating collapse 
by Shibata & Shapiro (2002) for n = 3 stars and 
by Duez et al. (2004) for n = 1 stars. In both cases 
the capture of fluid in orbits with j > jisco was 
demonstrated to be negligible. This criterion sug- 
gests a simple iterative scheme for calculating the 
final mass and spin of the hole and disk from the 
initial stellar density and angular momentum pro- 
file. First, guess the mass and spin of the hole, M/i 
and Jh- For our initial guess, we shall take a black 
hole that has consumed all the mass and angular 
momentum of the star, so that M^/M = 1 and 
Jh/Ml = {J/M'^)„if Next, use the initial stellar 
density and angular momentum profiles to "cor- 



rect" this guess by calculating the escaping mass 
and angular momentum of the outermost envelope 
with specific angular momentum exceeding jisco • 
We note that the value jisco depends on M/, and 
Jfi- We then "correct" the black hole mass and 
spin by deducting the values of the escaping mass 
and angular momentum of the envelope material 
from the guessed values of Mh and Jh- We recom- 
pute jisco for the "corrected" black hole mass and 
spin, and repeat the calculation of the escaping 
envelope mass and angular momentum until con- 
vergence is achieved. The calculation described 
exploits the theorem that for an axisymmetric dy- 
namical system, the specific angiilar momentum 
spectrum, i.e., the integrated baryon rest- mass of 
all fluid elements with specific angular momentum 
j less than a given value (e.g., iisco) is strictly 
conserved in the absence of viscosity (Stark & Pi- 
ran 1987). Any viscosity, if present, is expected 
to be unimportant on dynamical timescales, as re- 
quired by the theorem. 

For a Kerr black hole of mass Mh and spin pa- 
rameter a = Jh/Mh, the value of j'lsco is given 
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by 



Jisco = 



(20) 



where rms is the ISCO given by 
rms = M^S + ^2 - [(3 - Zi)(3 + Zi + 2Z2)]i/2}, 

where 



(21) 



a 



1/3 



1/3 



-) 



and 



1/3- 



1/2 



+ 



(22) 



(23) 



(see, e.g., Shapiro & Teukolsky 1983). Clearly, the 
infalling gas corotates with the black hole. 

The mass of the escaping matter in the envelope 
with j > jisco is given by 



AM : 



// 



2TTwdwp 



(24) 



where the density is given by cq. (13) and the 
quadrature is performed over all cylindrical shells 
in the star with cylindrical radii vo > ti^isco = 
(jisco/f^)^/^- Here we set Vt = flshcdd given by 
eq. (10). (The quantity computed in eq. (24) is 
actually the escaping rest-mass, but when the en- 
velope is nearly Newtonian, as we assume here, we 
can neglect the small difference between rest-mass 
and total mass-energy.) Defining w = w / Rp and 
z = z/Rp gives the nondimensional integral 



AM/M = 

X / / wdwdz 



1 



(25) 



(^2+^2)1/2 



Similarly, the angular momentum carried off by 
the escaping matter in the envelope is given by 



A J 



// 



(26) 



or 



Aj/M2 = ^f-") (^?i^'(a)i)^""'^ 



J J w^dwdz 



3/2 



Rp 
M 



(w2 + ^2)1/2 



1/2 



(27) 



where Rp/M is given in Table 2 and where once 
again the integral is performed over all cylindrical 
shells in the star with zu > zuisco = (jisco/^^)^^^- 
Eqns. (25) and (27) agree with Shapiro & Shibata 
(2002), eqns. (22) and (24), for n = 3. 

The mass and angular momentum of the black 
hole can then be determined from eqs. (25) and 
(27) according to 



Mh/M = l- AM/M 



and 



{J/M^ - AJ/M2 
(1 - AM/Mf 



(28) 



(29) 



Once the iteration of eqs. (20) - (23) with 
eqs. (25), (28) and (29) converges, the mass 
of the ambient disk can be found from 



Mdisk/M = AM/M. 



(30) 



Typically, convergence to better than 1% is 

achieved after only four iterations. 

The calculated parameters for the final black 
hole and disk are given in the last three columns 
of Table 2 for different values of n. The values for 
n = 3 are the most reliable, since the calculation of 
the density and angular momentum distributions 
in the outermost layers by means of a Newtonian 
Roche model is most accurate for such a soft equa- 
tion of state. For n = 3 we are able to compare 
the semi-analytic calculations with the numerical 
collapse simulation of Shibata & Shapiro (2002). 
Allowing for numerical error inherent in the simu- 
lation (< few percent), the agreement is quite good 
for the final masses of the black hole and disk, and 
reasonable, but less accurate, for the final black 
hole spin. The numerical simulation (third row in 
Table 2) used an initial model obtained by a nu- 
merical integration of the full general relativistic, 
equilibrium equations for the critical configuration 
at the onset of collapse. When the parameters for 
this model are also used in the analytic deriva- 
tion of the final black hole and disk (second row 
in Table 2) in place of the values determined by 
the variational principle (first row in Table 2), the 
agreement with the spin is much improved. 

The dependence of the black hole spin and disk 
mass on the polytropic index is summarized in Fig- 
ures 1 and 2. The spin parameter increases mod- 
erately with increasing n and increasing central 
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Fig. 1. — Mass- fraction of the disk formed fol- 
lowing collapse of a marginally unstable, supra- 
massive polytropc at the mass-shedding limit ver- 
sus polytropic index. The solid line uses criti- 
cal star parameters given by the variational treat- 
ment. The solid dots use critical parameters de- 
termined from numerical models of marginally un- 
stable, relativistic equilibrium stars at the mass- 
shedding limit (see Table 2). The solid triangle 
shows the value determined by a numerical simu- 
lation by Shibata & Shapiro (2002). 



concentration in the critical star. The disk mass 
is far more sensitive to the equation of state and 
increases very rapidly with increasing n. There 
is hardly any mass in the disk unless n is very 
close to 3, a result which is consistent with the 
dynamical simulations of Shibata, Baumgarte & 
Shapiro (2000) and Shibata (2003) covering the 
range 2/3 < n < 2. 

4. Application: Very Massive Star Col- 
lapse 

A very massive, uniformly rotating star {M/Mq ^ 
100) is characterized by the following properties: 

1. It is dominated by thermal radiation pres- 
sure; 



2. It is fully convective; 



Fig. 2. — Spin-parameter of the black hole formed 
following collapse of a marginally unstable, supra- 
massive polytropc at the mass-shedding limit ver- 
sus polytropic index. The curves and points are 
labelled as in Fig 1. 

3. It is governed by nearly Newtonian gravita- 
tion; 

4. It is described by the Roche model in the 
outer envelope. 

Such stars behave as polytropes with adiabatic 
index F = 1 -|- 1/n given by 



4/3 + /3/6 + 0(/32), 



(31) 



where Pr = is the radiation pressure and 

Pg = nkT/n is the gas pressure (see, e.g., Shapiro 
& Teukolsky 1983, Chapter 17). The ratio P <C 
1 is related to the radiation entropy per baryon 
Sr (« constant throughout the star) according to 
(3 — (4/yu)/(sr/fc), where /j, is the mean molecular 
weight, while the mass is related to the entropy 
according to 



Sr/k w 0.942(M/Mq). 



(32) 



Combining the above equations and evaluating 
/X « 2/(1 + 3X + 0.5Y) for a zero-metallicity, pri- 
mordial composition {X = 0.75 and Y = 0.25; Cy- 
burt. Fields and Olive 2003) likely to characterize 



8 



massive Population III objects, yields, to lowest 
order, 

^/^- " (3^- ^''^ 

Baumgartc & Shapiro (1999) have shown that, 
following cooling and contraction, such a mas- 
sive star will likely settle into rigid rotation and 
evolve to the mass shedding limit, assuming the 
viscous or magnetic braking timescale for angu- 
lar momentum transfer is shorter than the evolu- 
tion timescale in such an object (Bisnovatyi-Kogan 
et al. 1967; Zel'dovich & Novikov 1971; Shapiro 
2000; but see New & Shapiro 2001 for an alter- 
native scenario). Moreover, Newtonian simula- 
tions suggest that zero-metallicity rotating stars 
with M/Mq > 260 do not encounter the pair- 
production instability, so they collapse to black 
holes without exploding (Fryer, Woosley & Heger 
2001). The relation derived above between the 
mass M and polytropic index n for very massive 
stars then permits us to apply the results of the 
previous section to determine the final fate of these 
objects after they reach the onset of radial insta- 
bility and collapse. For example, if the star has 
a mass M/Mq > 10^, then 3 - n < 0.1. So ac- 
cording to Table 2 and Figs 1 and 2, the spin of 
the resulting hole will be moderate, J/j/M^ ^0-5, 
and the mass of the disk leftover from the collapse 
will be nonneghgible, Md/M > lO"^. The rela- 
tivistic simulations of Shibata & Shapiro 2002 for 
pure n = 3 polytropes serve to confirm these pre- 
dictions in the limiting regime 3 — n <C 1. 

5. Summary and Conclusions 

We have employed a simple analysis to de- 
termine the effect of the stiffness of the equa- 
tion of state on the fate of the collapse of a 
marginally unstable, relativistic polytrope spin- 
ning uniformly at the mass-shedding limit. We 
have used a variational principle and an approxi- 
mate, post-Newtonian energy functional to deter- 
mine the key parameters defining the structure of 
the critical progenitor star. We compared these 
parameters with the results of detailed nummeri- 
cal model calculations of stationary configurations 
for select cases. We then employed a Roche model 
to obain analytic expressions for the density and 
angular momentum profiles in the envelope of the 
marginally stable star. We substituted these pro- 



files into quadratures that determine the fractions 
of the stellar mass and angular momentum which 
escape capture by the central black hole assumed 
to form during the collapse. The fraction of the 
progenitor mass and spin which go into the hole 
versus the ambient disk is then iterated until con- 
vergence is achieved. 

We find that for the stars treated here, the mass 
fraction in the disk is about 10% for an n = 3 poly- 
trope and decreases rapidly as n decreases and the 
equation of state stiffens. The results are in agree- 
ment with the numerical simulations in 3 -|- 1 by 
Shibata, Baumgarte & Shapiro (2000) for n = 1 
and simulations by Shibata (2003) in axisymmetry 
for 2/3 < n < 2, which show that the mass frac- 
tion in a disk, if present at all, is less than 0.1% of 
the total mass. For the special case of n = 3, the 
results are also in agreement with the simulation 
performed by Shibata & Shapiro (2002) in axisym- 
metry, as has been discussed previously by Shapiro 
and Shibata (2002). The spin parameter of the 
black hole is less sensitive to the stiffness of the 
equation of state, decreasing from J^/M^ « 0.75 
for n = 3 to Jh/M^ « 0.39 for n = 2.5. 

The approach outlined here is more reliable for 
soft equations of state. In this limit the marginally 
unstable configuration is nearly Newtonian and 
the post-Newtonian energy functional describing 
the bulk of the mass, the Newtonian Roche model 
describing the envelope and the quadrature de- 
termining the escaping mass-energy fraction (as- 
sumed to equal the escaping rest-mass fraction) 
all become better approximations. Interestingly, 
this is precisely the limit where fully relativistic 
numerical simulations become more taxing, due to 
the large dynamic range and very long time inte- 
grations required for a hydrodynamic calculation. 

While the semi-analytic approach presented 
here to predict the final black hole and disk masses 
and spins was applied to treat the collapse of 
stars with uniform rotation, the method can be 
used equally well to treat the collapse of unstable 
stars with differential rotation. Recently, Duez et 
al. (2004) performed fully relativistic numerical 
simulations of hypermassive stars with apprecia- 
ble differential rotation. Hypermassive stars may 
form from the merger of binary neutron stars or 
from rotating core collapse in supernovae (Baum- 
garte, Shapiro & Shibata 2000). Magnetic braking 
or viscous damping of differential rotation in such 
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stars can drive them unstable to collapse on secu- 
lar timescales. Duez et al. (2004) performed sim- 
ulations to demonstrate this effect in the case of 
viscosity and showed that the final black hole and 
disk parameters are in agreement with the values 
predicted by the method outlined in Section 3. In 
particular, they showed that for rapid differential 
rotation characteristic of hypermassive stars, the 

disk mass fraction is typically large (10% 20% 

of the total initial mass) for stiff equations of state 
with n = 1. We therefore conclude that, in gen- 
eral, the parameters of the final black hole and 
the disk formed during the collapse of an unstable 
star depend both on the equation of state and on 
the degree of differential rotation. 
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